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ABSTRACT 


The  likely  extent  of  retardation  exerted  by  solvent  friction  upon  the  rates 
of  activated  electron-transfer  (ET)  processes  is  explored  and  evaluated  with 
reference  to  some  extant  solvent-dependent  data  for  metallocene  ET  self 
exchanges.  Barrier-crossing  frequencies  extracted  from  the  experimental  kinetic 
and  barrier  data  for  suitably  adiabatic  reactions  are  compared  with  solvent 
inertial  frequencies  (i.e.  the  zero-friction  limit)  as  estimated  from  currently 
available  analytic  expressions.  The  additional  extent  of  rate  retardation  seen 
in  going  from  low-friction  to  ostensibly  strongly  overdamped  solvents  as  deduced 
in  this  manner  is  seen  to  be  substantially  muted  in  comparison  with  the 
predictions  of  conventional  Debye-continuum  approaches,  although  follows  a 
solvent  functionality  that  is  nevertheless  roughly  in  accord  with  the  latter. 
The  likely  importance  of  more  rapid  dynamics  associated  with  short-range 
(molecular)  solvation  is  addressed  in  the  light  of  these  findings. 
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INTRODUCTION 

Elucidating  the  various  roles  of  the  molecular  solvent  on  the  kinetics  of 

activated  rate  process  which  occur  in  such  media  is  a  quest  that  is  not  only  of 

broadbased  fundamental  interest  and  importance,  but  is  seemingly  of  sufficient 

nobility  to  have  attained  "Holy  Grail-like"  status  in  the  minds  of  many 

experimentalists  and  some  theoreticians.  Perhaps  inevitably,  earlier  treatments 

and  discussions  focussed  on  general  thermodynamic  and  activation  energetic 

aspects.^-  Beginning  in  the  late  1970' s,  and  swelling  (as  in  blood  clots)  during 

the  last  decade  up  to  the  present  time,  has  been  a  substantial  upsurge  of 

interest  in  the  influence  of  the  solvent  upon  reaction  dynamics .  Progress  has 

been  fueled  in  particular  by  extensive  developments  in  the  theoretical 

description  of  so-called  "solvent  friction,"  whereby  the  net  dynamics  for 

activated  barrier  crossing  (and,  more  generally,  along  the  reaction  coordinate) 

are  slowed  below  the  u«  ;al  values  prescribed  by  transition-state  theory  (TST) 
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treatments  by  dissipative  energy-transfer  effects  involving  the  solvent  "bath." 

The  precise  physical  meaning  and  interpretation  of  solvent-friction  effects 

depend  somewhat  on  the  type  of  chemical  process  being  considered,  specifically 

on  the  nature  of  the  reaction  coordinate.  Nonetheless,  the  usual  notion 

perceives  the  presence  of  an  impediment  to  progress  along  the  reaction  coordinate 

as  arising  from  irreversible  energy  dissipation  from  the  reactive  modes  to  the 

surrounding  solvent.  The  reaction-coordinate  motion  is  then  no  longer  smooth  and 

unidirectional  but  is  characterized  instead  by  "fits  and  starts",  commonly  termed 

"overdamped  motion",  perhaps  even  involving  several  recrossings  back  and  forth 

2 

within  the  barrier-top  region  before  the  reaction  is  finally  consummated. 

One  important,  although  historically  not  the  first,  class  of  activated 
reactions  for  which  the  notion  of  solvent-friction  effects  has  been  considered 
extensively  is  electron  and  related  charge  transfers.  (See  refs.  3-7  for  some 


recent  overviews).  Such  systems  involve  solvent  motion  as  a  central  component 
of  the  reaction  coordinate,  specifically  in  the  marked  changes  in  the 
polarization  of  the  solvent  dipoles  commonly  required  for  the  occurrence  of 
charge  transfer  between  reactive  sites.  This  collective  dielectric  polarization 
is  coupled  inherently  with  charge  transfer.  However,  the  rapidity  of  electronic 
compared  with  nuclear  motion  obliges  the  latter,  including  solvent 
reorganization,  to  occur  at  least  partly  prior  to  the  instant  of  charge  transfer 
itself,  thereby  giving  rise  to  the  nuclear  activation-energy  barrier.  In  cases 
where  solvent  reorganization  constitutes  the  predominant  contribution  to  the 
free-energy  barrier,  this  collective  solvent  motion  necessarily  determines  the 
dynamics  of  approaching,  and  possibly  also  surmounting,  the  activation  barrier. 
As  a  consequence,  one  can  anticipate  that  such  solvent  "dielectric  friction" 
could  exert  a  significant  or  even  substantial  influence  upon  the  net  barrier¬ 
crossing  dynamics,  and  hence  on  the  ET  reaction  rate. 

This  intriguing  expectation  was  indeed  initiated  and  then  fueled  by  a 

number  of  theoretical  contributions  appearing  in  the  early/mid  1980 's  (refs.  8-16 

are  a  very  incomplete  selection  of  key  papers)  .  The  emergence  of  this 

theoretical  framework  spawned  a  large  number  of  experimental  studies  aimed 

variously  at  probing  the  nature  of,  and  the  extent  to  which,  collective  solvent 

polarization  can  influence  the  dynamics  of  ET  and  related  processes.  Broadly 

speaking,  the  experimental  inquiries  fall  into  four  distinguishable  (although 

partly  overlapping)  categories.  The  first  class  (type  I)  involves  the  evaluation 

of  time-dependent  fluorescence  Stokes  shifts  (TDFS)  for  solute  chromophores 
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forming  suitable  charge-transfer  excited  states.  ’  In  optimal  cases,  these 
ultrafast  (>  50  ps)  laser-induced  measurements  can  probe  directly  the  real-time 
dynamics  of  polar  solvent  relaxation  around  »  newly  formed  solute  dipole.  The 
second  (type  II)  concerns  a  related  use  of  ultrafast  pulse  lasers  to  follow 
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intramolecular  ET  processes  emanating  from  photoexcited  states.  So  far, 

however,  almost  all  of  the  ET  reactions  studied  in  this  fashion  involve  only 
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small  free-energy  barriers,  AG*  5  kBT  (kB  is  Boltzmann's  constant).  ’ 

The  remaining  two  experimental  categories  constitute  attempts  to  extract 
information  concerning  dynamical  solvent  effects  for  activated  ET  processes  (say, 
for  AG*  >  5kBT) .  The  most  common  approach  (type  III)  involves  the  acquisition 
of  solvent-dependent  kinetic  data  for  bimolecular  outer-sphere  electron-exchange 
reactions,  especially  symmetrical  self-exchange,  and  for  related  reactions  at 
electrode  surfaces.**  The  kinetics  of  the  former  are  most  commonly  evaluated  by 
means  of  magnetic  resonance  (nmr,  epr)  line-broadening  methods.  These  reaction 
systems  have  the  obvious  virtues  of  featuring  significant  (commonly  5-10  kBT) 
solvent  reorganization  barriers,  thereby  being  truly  activated  processes,  and 
provide  opportunities  to  alter  systematically  the  nature  of  the  solute-solvent 
interactions,  reactant  electronic  structure  and  vibrational  barriers,  etc. 
Indeed,  a  substantial  amount  of  information  regarding  solvent  dynamical  effects 
can  in  principle  be  gleaned  from  such  kinetic  data.  Despite  the  clear 
fundamental  importance  of  activated  ET  processes,  however,  the  extraction  of  the 
required  information  on  barrier-crossing  dynamics  is  often  problematical ,  leading 
to  the  publication  of  a  number  of  rate  analyses  with  accompanying  (sometimes 
wide-ranging)  conclusions  that  the  connoisseur  (or  even  casual  observer)  might 
justifiably  regard  as  specious  at  best. 

The  core  of  the  problem  concerns  the  need  to  separate  the  dynamical 
(preexponential  factor)  and  energetic  (barrier  height,  AG*)  contributions  to  the 
measured  solvent-dependent  rate  constants.  In  most  cases,  the  extraction  of  the 
desired  dynamical  information  from  rate  data  obliges  one  to  resort  to  the  use  of 
theoretical  estimates  for  AG*,  usually  having  uncertain  quantitative  validity. 
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An  additional  uncertainty  in  these  conventional  analyses  is  that  the  geometry  as 
well  as  stability  of  the  precursor  complexes  for  bimolecular  (or  electrochemical) 
outer-sphere  processes  is  often  unknown. 

The  conclusions  regarding  the  role  of  solvent  dynamics  reached  in  this 
fashion  are  therefore  often  critically  dependent  on  the  validity  of  the 
(suspected  oversimplified)  theoretical  models  utilized  to  separate  the  dynamical 
and  energetic  factors.  These  type  III  tactics  have  consequently  tended  to  be 
regarded  as  a  "poor  cousin"  in  comparison  with  the  more  direct  (and 
sophisticated-looking)  insight  into  real-time  dynamics  obtained  from  approaches 
I  and  II.  Nevertheless,  in  some  cases  sufficient  information  on  barrier  height 
and  related  factors  can  be  obtained  so  to  allow  more  reliable,  and  even 
quantitative,  deductions  to  be  made  regarding  the  solvent-dependent  barrier¬ 
crossing  dynamics  for  bimolecular  ET  reactions  (vide  infra).** 

Most  of  the  uncertainties  faced  with  such  dynamical  analyses  for  activated 
ET  reactions  can  also  be  circumvented,  at  least  in  principle,  by  evaluating 
solvent-dependent  rate  data  for  intramolecular  (preferably  symmetric)  ET 
reactions  (labelled  here  type  IV).  The  evaluation  of  both  the  unimolecular  ET 
rates  and  the  barrier  heights  (the  latter  from  optical  ET  energies)  would  provide 
the  preferred  route  to  the  separation  of  dynamical  and  energetic  factors. 
Unfortunately,  however,  relatively  few  type  IV  systems  have  been  scrutinized 
usefully  so  far  in  the  context  of  solvent  dynamics. ^ 

These  overall  issues  have  been  reviewed  in  some  detail  recently  by  the 

ga 

present  author.  Broadly  speaking,  the  extant  experimental  information  lends 
support  to  the  earlier  theoretical  predictions  that  overdamped  solvent  motion 
("solvent  friction")  may  influence  and  even  control  the  ET  barrier-crossing 
dynamics,  as  prescribed  most  simply  by  the  longitudinal  relaxation  time,  rL,  in 
a  number  of  cases.  However,  several  interesting  and  even  crucial  issues  remain 
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largely  undiscussed  as  well  as  unresolved.  One  nagging  question  involves  the 
differences  between  the  dynamics  of  short-range  (molecular-scale)  solvation, 
which  should  provide  a  dominant  part  of  the  net  solvent  friction,  and  the 
inherently  long-range  solvent  dynamics  considered  in  the  dielectric-continuum 
models  which  are  commonly  utilized  in  data  analyses.  While  the  latter  (described 
for  Debye  solvents  simply  by  r^*1)  correlates  with  the  solvent-dependent  barrier¬ 
crossing  rates  in  a  number  of  experimental  systems,^  one  has  good  reason  to  be 
suspicious  of  the  common  inference  that  rL-1  "describes  the  relevant  dynamics." 
A  related  problematic  issue  is  the  widespread  observation  (or  claims!)  of  such 
rate-rL_1  correlations,  especially  for  electrochemical  reactions, even  though 
established  electron-tunneling  treatments  predict  their  occurrence  only  as 
limiting  cases  for  ET  reactions  displaying  unusually  strong  donor— acceptor 
electronic  coupling. ^ 

The  purpose  of  the  present  paper  is  to  address  in  straightforward 
descriptive  fashion  these  questions  concerning  ET  barrier-crossing  rates  in  the 
light  of  current  notions  concerning  molecular  solvent  dynamics,  that  is, 
concerned  especially  with  short-range  solvation.  The  basic  tenet  involves  the 
use  of  selected  experimentally  derived  data  for  barrier-crossing  frequencies  in 
comparison  with  expectations  from  theoretical  models  in  order  to  assess  the 
validity  of  conventional  descriptions  of  solvent-friction  effects  in  ET 
reactions . 

Extraction  of  Solvation  Dynamical  Information  for  Activated  ET  Processes 

We  consider,  as  before,  homogeneous-phase  ET  reactions  having  the  general 

form 


0xx  +  Red2  »  Redj  +  0x2 


(1) 
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where  Ox  and  Red  refer  to  oxidized  and  reduced  forms  of  a  given  redox  couple. 
We  are  concerned  primarily  with  self-exchange  reactions,  where  Oxj  -  0x2  and  Red2 
-  Redlt  so  that  the  free-energy  driving  force,  AG",  is  zero,  and  the  reaction  is 
chemically  as  well  as  energetically  symmetrical.  The  Ox  and  Red  sites  may  reside 
on  separate  solute  species  (intermolecular  ET) ,  or  may  be  bound  chemically  to 
each  other  (intramolecular  ET) .  Experimental  kinetic  results  for  only  the  former 
type  (III  above)  will  be  considered  here. 

The  role  of  the  dynamics  of  nuclear  modes,  including  those  associated  with 
solvent  reorganization,  upon  the  rate  constant  for  electron  exchange,  k„,  is 
commonly  expressed  in  terms  of  the  nuclear  barrier-crossing  frequency,  (s_1), 
by  means  of  the  simplified  expressions^ 

(2a) 

where 

Kt  “  *ei  ‘'n  exp  (-AGj/kBT)  (2b) 

Here  AGj  is  the  intrinsic  free-energy  barrier,  *el  is  the  electronic  transmission 
coefficient,  and  Kp  is  the  prequilibrium  constant  describing  the  statistical 
probability  of  finding  the  reactant  pair  in  an  internuclear  configuration 
suitable  for  electron  transfer.  For  intramolecular  ET,  the  desired 
"unimolecular"  rate  constant  k,t  (s_1)  is  usually  measured  directly,  whereas  for 
bimolecular  reactions  it  is  necessary  to  estimate  k,t  from  the  measured  rate 
constant  k#x  by  means  of  formal  expressions  such  as  Eq.  (2a). 

A  central  fundamental  issue  that  has  incurred  much  discussion  is  the  manner 
and  extent  to  which  overdamped  solvation  dynamics  ("solvent  friction")  may  limit 
the  measured  rate  constant  by  affecting  An  important  aspect  of  this 
question  concerns  the  relative  roles  of  donor-acceptor  electronic  coupling  and 
nuclear  dynamics  in  determining  the  overall  preexponential  factor  A.t  (-  x#1  wn) . 
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Simply,  for  systems  where  tnc  electronic  coupling  (as  prescribed  by  the 

electronic  matrix  coupling  H12)  is  small  (typically,  say,  <  0.5  kJ  mol-1’),  A,t 

will  be  dominated  by  the  electron-tunneling  probability  within  the  transition 

state,  and  the  nuclear  dynamics  (i/n)  should  exert  little  or  no  influence  upon  A,t 

( "nonadiabatic  ET").  Only  for  systems  where  the  electronic  coupling  is 

sufficient  to  yield  adiabatic  pathways  (k#1  -*■  1)  will  nuclear  dynamics  constitute 

the  dominant  contribution  to  A.t.  The  latter  situation  corresponds  to  the 

reaction  dynamics  being  describable  in  terms  of  motion  along  a  single 

hypersurface  continuously  from  equilibrium  reactants  to  products.  Most 

significantly  from  an  experimental  perspective,  therefore,  the  crossover  from 

largely  nonadiabatic  to  adiabatic  behavior  is  expected  to  be  signaled  by  a  marked 

alteration  in  the  solvent-dependent  kinetics  in  that  the  reaction  dynamics  should 

18 

be  sensitive  to  the  solvation  dynamics  only  in  the  latter  case. 

In  our  laboratory,  this  prediction  turned  out  to  have  a  major  beneficial 

influence  on  our  efforts  to  disentangle  solvent  dynamical  effects  from  the  other 

factors  which  may  influence  the  rates  of  type  III  ET  processes.  Rate  data  for 

a  series  of  metallocenium-metallocene  (Cp2M+/0,  where  Cp  -  cyclopentadienyl  and 

derivatives)  self-exchange  reaction  were  obtained  in  a  number  of  solvents  with 

19-22 

suitably  varied  dynamical  properties.  [The  reasoning  behind  the  choice  of 

these  systems  as  "model  probe  reactants"  for  solvent  dynamical  studies  is 

detailed  in  the  original  papers.  Nonetheless,  worthy  of  mention  here  is  the 

availability  of  experimental  estimates  of  AG*  from  optical  ET  energies  for  mixed- 

24  25 

valence  Cp2Fe+/0  dimers,  ’  which  facilitate  the  (hopefully)  reliable  separation 

19 

of  energetic  and  dynamical  factors  (vide  supra).]  While  not  evident  initially, 
subsequent  work  Identified  clearly  a  crossover  in  the  solvent-dependent  kinetics 
from  the  relatively  slow  reactions  (exemplified  by  Cp2Fe'f/0)  ,  which  display 
virtually  solvent  -  independent  dynamics,  to  more  facile  processes  (such  as  for 
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Cp2Co+/0  derivatives)  that  yield  approximately  rL_1-dependent  dynamics. 

Moreover,  this  transitional  behavior  was  found  to  be  consistent  with  theoretical 

22  23 

H12  estimates  for  Cp2M+/0.  ’ 

Besides  the  intrinsic  significance  of  the  observed  nonadiabatic— adiabatic 

crossover,  the  results  provide  an  unusual  (even,  we  believe,  so  far  unique) 

opportunity  to  assess  the  validity  of  the  kinetic  data  analysis  utilized  for 

extracting  solvent-dependent  dynamical  information.  Most  simply,  one  can 

crosscheck  the  absolute  magnitude  as  well  as  the  solvent  dependence  of  the  A*t 

values  obtained  from  the  experimental  kax  and  AG*  values  by  using  Eq .  (2)  or 
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related,  more  sophisticated,  expressions.  As  a  consequence,  we  have  some 

confidence  (so  far!)  in  at  least  the  semiquantitative  validity  of  estimates  of 
A.t  and  related  quantities  extracted  from  Cp^'0  self-exchange  data.  We  cannot 
make  a  similar  claim  for  any  other  type  III  system  studied  to  date. 

Estimation  of  Inertial  Solvation  Dynamics 

Given  this  scenario,  It  should  be  instructive  to  compare  further  this 
experimentally  derived  information  on  solvent-dependent  barrier-crossing 
dynamics,  published  now  four  years  ago,  with  the  expectations  of  current, 
especially  molecular-based,  models  of  solvation  dynamics.  Of  particular 
significance  in  this  regard  is  the  extent  to  which  overdamped,  rather  than 
inertial,  solvent  motion  controls  the  barrier-crossing  dynamics  in  adiabatic  ET 
processes.  The  latter  situation  corresponds  to  the  occurrence  of  the  transition- 
state  theory  (TST)  limit,  where  solvent-friction  effects  are  absent.  In  this 
case,  t/D  will  be  determined  instead  by  the  "free"  rotational  dynamics  of 
individual  solvent  dipoles,  so  that  vn  -  w0/2*,  where  u0  is  this  inertially 
limited  frequency. 

An  interesting  means  of  gauging  the  importance  of  solvent-friction  effects 
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is  to  assess  the  degree  to  which  the  experimental  solvent-dependent  values  of 

<c#1vn  for  adiabatic  ET  processes  (where  x#1  -*  1)  fall  below  the  anticipated 

solvent  inertial  frequencies.  To  be  numerically  worthwhile,  such  a  comparison 

would  desire  the  receipt  of  relatively  reliable  estimates  of  u0/2n  (say,  within 

9  ii 

twofold).  Following  earlier  treatments,  '  a  pair  of  analytic  expressions  for 
26 

w0  can  be  derived.  The  first  can  be  expressed  as 

u0  -  uf[(2 e,  +  e«)/3e«gk]  (3) 

Here  es  and  ea  are  the  static  and  "infinite-frequency"  dielectric  constants,  gk 
is  the  Kirkwood  g-f actor,  and  u£  is  the  "free  rotor"  frequency  given  by 
(2kBT/I)%,  where  I  is  the  solvent  molecular  moment  of  inertia  (as  extracted  from 
gas-phase  microwave  spectra).  The  square-bracketed  term  in  Eq.  (3)  describes  how 
the  effective  solvent  rotational  frequency  is  accelerated  by  the  presence  of  the 
surrounding  dielectric  medium  (vide  infra) . 

The  second  expression  for  w0  can  be  obtained  from  Eq .  (72)  of  ref.  9,  which 

,  .  26 
can  be  written  as 

w0  -  wr(3y <,/(<,  -  1)]'  (A) 

The  reduced  concentration  of  solvent  dipoles,  y,  is  given  by 

y  -  4wpp2/9kBT  (4a) 


where  p  is  the  solvent  number  density  and  p  is  the  dipole  moment.  The  close 

relationship  between  Eq.  (4)  and  the  alternative  expression  for  o>0,  Eq.  (3),  is 

27 

exposed  by  writing  the  Kirkwood  g-factor,  appearing  in  the  former  relation,  as 


„  .  iiSbl  !&.=&,)  C Zi,±-jul 
Bk  4ir p2  <s(c0  -i-  2)2 


(5) 


Equations  (3)  and  (4)  are  equivalent  when  gk  is  given  by  Eq.  (5),  in  the  limit 
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where  1.  This  circumstance  is  tantamount  to  neglecting  the  electronic 

polarization  of  the  solvent  medium. 

An  illuminating  discussion  of  the  factors  influencing  solvent  inertial 

dynamics,  utilizing  essentially  Eq.  (4),  has  been  presented  recently  by 
2  8 

Maroncelli  et  al .  They  employed  (and  derived)  Eq.  (4)  expressed  in  the  form 

«o  “  a  (6) 

where  the  coefficient  a  is  given  by 

q  -  (47rp/j2/3kBT)  (1  -  es-1)-1  (6a) 

Given  that  es  ►  1  for  polar  solvents,  so  that  (1  -  «B_1)  »  1,  this  formulation 

emphasizes  how  the  effect  of  the  dielectric  environment  on  the  inertial  dynamics 

is  predicted  to  be  dominated  simply  by  the  solvent  dipole  density  as  gauged  from 

the  ppz  term.  Maroncelli  et  al  found  support  for  the  validity  of  Eq.  (6a)  from 

the  approximate  agreement  between  the  w0  values  predicted  in  this  fashion  and 

those  obtained  from  MD  simulations  for  several  solvents,  including  acetonitrile 
28 

and  water.  Since  the  latter  calculations  refer  largely  to  solvation  dynamics 

involving  spherical  solutes,  they  would  appear  to  have  direct  relevance  for 

predicting  the  inertial  dynamics  associated  with  short-range  solvation  in 

addition  to  the  "bulk-phase"  solvent  to  which  s<’ch  dielectric  continuum  models 

usually  refer.  One  therefore  has  a  firm  basis  for  confidence  in  the  usefulness 

of  Eq.  (4)  for  estimating  inertial  solvation  dynamics  in  activated  ET  processes. 

At  first  glance,  a  limitation  of  Eqs .  (4)  and  (6a)  would  appear  to  be  the 
neglect  of  the  high  (optical)  frequency  dielectric  constant,  which  is  included 
in  the  alternative  formulation  Eq .  (3).  Combining  Eqs.  (3),  (5)  and  (6)  yields 

€,(€.+2)3 

9e„(e,-0 


a=  (4rcpp2/3  kBT) 


(6b) 
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where  the  term  in  square  brackets  now  replaces  the  (1  -  t,"1)"1  term  in  Eq .  (6a). 
Interestingly,  however,  the  former  term  is  consistently  close  to  unity  (within 
ca  20%)  for  values  of  e4  and  (>  20  and  2-4,  respectively)  typical  for  polar 
solvents,  so  that  the  inclusion  of  the  electronic  polarizability  of  the  solvent 
molecules  (via  <*)  exerts  little  numerical  influence  upon  the  estimation  of  w0. 
The  central  role  of  the  dipole  density  in  prescribing  the  effect  of  the 
surrounding  dielectric  medium  upon  the  inertial  frequency  is  therefore  an 
intriguing  hallmark  of  both  formulations,  Eqs .  (3)  and  (4). 

Assessment  of  Rate  Retardation  by  Solvent  Friction 

On  the  basis  of  the  foregoing,  the  comparison  between  the  solvent-dependent 
barrier-crossing  frequencies,  n,1i/n,  derived  from  the  Cp2M+/0  rate  and  energy 
barrier  data  with  inertial  frequencies  estimated  from  Eqs.  (3)  and  (4)  should 
provide  at  least  semiquantitative  insight  into  the  degree  to  which  overdamped 
dynamics  (solvent  friction)  may  diminish  the  rate  below  the  adiabatic  TST  limit. 
To  this  end,  Table  I  displays  Kali/n  values  evaluated  for  three  metallocene  self¬ 
exchange  reactions  in  eleven  solvents,  listed  alongside  the  corresponding  solvent 
inertial  frequencies,  w0/2tt,  as  estimated  from  Eq .  (3).  (See  refs.  20  and  21  and 
footnotes  to  Table  I  for  data  and  kinetic  analysis  used  to  obtain  values; 

see  footnotes  to  Table  II  and  Fig.  4  of  ref.  29  for  w0  calculational  details.) 
The  three  metallocene  couples  consist  of  a  pair  of  cobaltocene  systems,  Cp2Co+/0 
and  CpfCo+/0  (Cp*  -  pentamethylcyclopentadienyl ,  Cp*  -  carboxymethyl  derivative), 

that  follow  largely  adiabatic  pathways,  and  one  example,  CpzFe+/0,  that  displays 

22 

the  hallmarks  of  reaction  nonadiabaticity .  The  latter  behavior  is  evidenced 
most  clearly  from  the  observed  insensitivity  of  the  k.iv„  values  for  the 
ferrocene  couple  to  the  solvent  dynamics,  as  assessed  conventionally  by  the  rL~l 
values  that  are  also  listed  in  Table  I.  In  contrast,  the  cobaltocene  reactions 
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display  markedly  (ca  10-100  fold)  larger  barrier-crossing  frequencies;  moreover, 

the  values  correlate  roughly  with  the  anticipated  solvent  dynamics,  at 

least  in  the  "Debye"  (.first  nine)  solvents  listed  in  Table  I.  This  correlation 

is  in  approximate  concordance  with  the  expectation  that  solvent  dynamics  controls 

the  adiabatic  barrier-crossing  frequency,  since  according  to  the  simplest 

overdamped  dielectric-continuum  models  one  anticipates  that^  uD  -  rL-1  (but  see 

below).  We  have  discussed  at  length  elsewhere  this  issue  with  respect  to  the 

21  22 

metallocene  reactions.  ’  In  particular,  a  more  detailed  analysis  of  the 

variations  in  the  rate-solvent  dependencies  for  the  different  metallocene  self- 

22 

exchanges  enabled  estimates  of  the  electronic  coupling  to  be  deduced.  Thus  for 
the  largely  adiabatic  cobaltocene  couples,  H12  -  0.2  to  5  kJ  mol-1,  whereas  for 
the  more  nonadiabatic  ferrocene  systems,  H12  -  0.4  kJ  mol-1. 

A  close  comparison  of  the  solvent-dependent  cobaltocene  values  with 

the  corresponding  w0/2ir  as  well  as  rL-1  values  in  Table  I,  however,  suggests  that 
a  significantly  different  picture  of  the  role  of  solvent  dynamics  may  be  more 
appropriate.  Specifically,  the  values  for  both  cobaltocene  reactions  lie 

almost  uniformly  within  ca  twofold  of  the  predicted  inertial  barrier-crossing 
frequencies,  u0/2n ,  in  each  solvent,  implying  that  the  actual  frictional 
retarding  effect  is  no  more  than  small  or  moderate.  This  observation,  of  course, 
must  be  tempered  by  the  recognition  of  the  likely  uncertainties  (probably  around 
twofold  or  so)  in  estimating  as  well  as  w0.  Nevertheless,  when  fueled  by 

the  likely  reliaV  lity  of  the  inertial  frequency  estimates  as  noted  above,  one 
seemingly  has  grounds  for  doubting  the  oft-perceived  dominance  of  overdamped, 
rather  than  inertial,  solvent  dynamics  in  adiabatic  ET  barrier  crossing  even  for 
systems  such  as  the  metallocenes  where  solvent  reorganization  forms  the 
predominant  barrier  component.  Thus  taking  the  x#1i/n  and  u0/2n  estimates  at  face 
value,  the  degree  to  which  solvent  friction  retards  the  barrier-crossing 
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frequency  below  the  inertial  (TST)  limit  would  appear  to  be  no  more  than  2-3 
fold,  even  in  supposed  "slow  well-overdamped’’  solvents  such  as  benzonitrile ,  TMU, 
and  HMPA. 

It  is  appropriate,  however,  to  review  briefly  the  basis  of  our  earlier 

suggestion  that  overdamped  dynamics  provide  a  predominant  influence  upon  for 
6  19_22 

such  ET  reactions.  ’  The  values  obtained  for  the  adiabatic  (or  near- 

adiabatic)  reactions  decrease  by  10-20  fold  in  going  from  dynamically  fast  (and 

low  friction)  solvents  such  as  acetonitrile  and  water  to  "slow  overdamped"  media 

such  as  TMU  and  HMPA  (Table  I).  While  the  corresponding  variations  in  the  rL_1 

values  are  somewhat  larger,  20-30  fold,  such  sub-unit  log  -  log  rL_1  slopes 

are  expected  even  on  the  basis  of  simple  overdamped  solvent  dynamical  models  on 

several  grounds. Most  generally,  the  estimation  of  solvent-dependent 

values  presuming  a  fixed  precursor  stability  (Kp) ,  as  was  done  in  assembling 

Table  I,  is  an  oversimplification.  Kp  is  actually  anticipated  to  diminish 

progressively  as  i/n  increases  as  a  result  of  systematic  alterations  in  the 

distribution  of  encounter  geometries  co:lt  ibuting  to  the  measured  bimolecular 

rate  constant. (This  effect  can  partly  be  accounted  for  within  the 

framework  of  the  "fixed  Kp"  analysis  by  noting  that  k81  for  near-adiabatic 

18 

reactions  should  decrease  as  vn  increases.  )  A  formally  more  precise  analysis 

involves  the  comparison  of  observed  rate-solvent  dependencies  with  the 

corresponding  calculated  behavior  that  accounts  for  precursor  geometry  and  other 

complicating  features:  such  a  procedure  applied  to  the  metallocene  systems 

6  22 

yielded  valuable  insight  into  the  role  of  electronic  coupling.  ’  In  addition, 

the  presence  of  reactant  vibrational  contributions  to  the  ET  barrier  can  also 

12 

diminish  the  x8li/n  -  rL-1  dependence,  even  though  the  observed  attenuation 

30 

appears  typically  to  be  milder  than  predicted. 

Several  theoretical  expressions  emanating  from  the  dielectric-continuum 
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model  provide  estimates  of  the  extent  of  solvent  frictional  retardation  of  the 

ET  barrier-crossing  dynamics  for  certain  specific  forms  of  the  free-energy 

barrier.  The  first  such  description  of  the  interplay  between  overdamped  and 

9 

inertial  solvent  dynamics  was  given  by  Calef  and  Wolynes.  Their  expression  for 

the  adiabatic  barrier-crossing  frequency  for  electron-exchange  processes  can  be 
34 

written  as 

vfl=  (2ti cr.) _1  [0 . 5  ♦  0 . 5  ( 1  — ~ —  •— ) Vfc] (7) 

cae,Y  Tj.2 

Here  7  -  (3y/(es-l)]  {Eq.  (4a)],  rrot  -  (Wf)"1,  and  c  is  a  so-called  barrier  "shape 
9 

factor."  For  cusp-like  barriers,  i.e.  H1Z  <  kBT,  this  quantity  is  given 

approximately  by  c  -  (kBT)/7rAG*)' .  In  the  limit  where  rrot  <  rL,  Eq.  (7)  reduces 

g 

to  the  well-known  expression  for  overdamped  dynamics  within  a  cusp  barrier 

vn  -  rL"1(AG*/4irkBT)'1  (8) 

In  the  inertial  limit,  i.e.  where  rrot  ►  rL,  we  recover  the  TST  formula 

‘'n  -  wo/2*  (9) 

The  "shape  factor"  c  in  Eq.  (7)  accounts  for  the  effects  of  barrier-top 

roundedness  brought  about  by  donor-acceptor  electronic  coupling.  Since  c 

increases  towards  larger  H12,  Eq.  (7)  predicts  that  the  barrier-crossing 

frequency  in  the  presence  of  solvent  friction  will  diminish  as  the  barrier  top 

35 

becomes  more  rounded  (i.e.  less  cusp-like).  This  barrier-shape  dependence  of 
the  reaction  dynamics,  which  is  absent  in  the  TST  (inertial)  limit,  arises 
physically  from  the  additional  "recrossings”  associated  with  overdamped  solvent 
motion  that  are  anticipated  as  the  barrier  top  becomes  broader. 

An  alternative  combined  formulation  of  overdamped  and  inertial  motion, 
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which  emphasizes  more  explicitly  the  influence  of  barrier-top  dynamics  on  the 
barrier-crossing  frequency,  has  been  described  by  Hynes. ^  This  treatment 
pictures  the  overall  adiabatic  rate  constant,  kA,  as  a  composite  of  the  rate  for 

crossing  the  barrier  top,  kb,  and  for  motion  within  the  parabolic  wells,  k*  , 

..  .  11 
according  to 

kA_1  -  V1  +  0.5  k/1  (10) 

While  solvent  friction  is  included  in  the  analytic  expressions  for  both  kb  and 

k*,  for  reasons  of  consistency  solvent  inertia  is  considered  only  in  the  former 

component.  Nonetheless,  for  Debye  solvents  the  Hynes  formulation  [Eq.  (10)] 

yields  comparable  numerical  results  to  the  Calef-Wolynes  expression  [Eq.  (7)]. 

A  key  advantage  of  the  former,  however,  is  that  it  is  designed  to  be  applicable 

in  the  presence  of  so-called  "frequency— dependent  solvent  friction. This 

circumstance  includes  not  only  that  encountered  necessarily  in  non-Debye 
29 

solvents  but  also  for  nonexponential  decay  dynamics  arising  from  molecular 

36 

solvation  effects.  Additionally,  the  Hynes  formulation  provides  a  more 

generally  useful  treatment  of  barrier-shape  effects  upon  the  barrier-crossing 

dynamics.  Thus  for  cusp  barriers,  kb  ►  k*  so  that  well  relaxation  dominates;  kb 

will  decrease,  however,  as  the  barrier  top  becomes  broader  so  that  eventually  kb 

5  k*,  whereupon  the  net  rate  will  diminish  to  the  point  where  it  is  controlled 

29 

by  solvent  dynamics  within  the  barrier-top  region. 

At  least  in  Debye  media,  both  formulations  predict  that  the  inertial 

dynamics  exert  an  important  limiting  influence  on  only  in  "low-friction"  media 

such  as  acetonitrile  and  acetone,  the  barrier-crossing  frequency  in  higher- 

friction  media  such  as  benzonitrile  [where  (w0/2x)  >  2rL-1]  being  controlled 

29  34 

almost  entirely  by  overdamped  solvent  dynamics.  ’  Given  that  the  near- 
adiabatic  cobaltocene-based  reactions  in  Table  I  necessarily  feature  substantial 
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barrier-top  roundedness  [H12  -  2-4  kJ  mol'1  (vide  supra)  ],  it  is  of  particular 

interest  to  consider  the  solvent-dependent  barrier-crossing  frequencies  predicted 

from  Eq.  (10)  under  these  conditions.  Listed  in  Table  1  are  values,  labelled 

y®,  calculated  from  the  Hynes  formulation,  Eq.  (10),  by  means  of  the  dielectric- 

continuum  model  using  the  u0  and  rL  values  listed  alongside,  with  kb  estimated 

for  H12  -  2.1  kJ  mol-1.  (This  matrix  coupling  value  is  appropriate  for  the 

22 

cobaltocene  reactions  considered  in  Table  I.)  For  the  non-Debye  solvents 
included  in  Table  I,  propylene  carbonate  (PC)  and  methanol,  the  distribution  of 
relaxation  times  indicated  in  the  dielectric  loss  spectra  is  taken  into  account 
when  estimating  kb  and  lq,.  Further  details  of  the  calculational  procedures  are 
given  in  ref.  29.  The  v®  values  listed  in  Table  I  are  markedly  (3-4  fold) 
smaller  than  r£l,  reflecting  the  predicted  sluggishness  of  surmounting  the 
barrier  with  a  well-rounded  top  by  means  of  overdamped  solvent  motion. 

Comparison  of  these  i/®  estimates  with  the  experimental  values  for  the 

cobaltocene  reactions  in  Table  I  shows  that  the  latter  are  typically  3-10  fold 
larger.  Recalling  our  observation  above  that  the  values  tend  to  lie  only 

twofold  or  so  below  the  corresponding  TST  limit,  w0/2»r,  one  may  conclude  that  the 
extent  of  solvent  friction  as  predicted  from  the  most  relevant  dielectric- 
continuum  approach  is  rather  more  pronounced  than  occurring  experimentally,  at 
least  for  the  metallocene  systems.  Stated  more  precisely,  the  experimental 
extent  of  rate  retardation  by  solvent  friction,  defined  as  kTST/k,  would  appear 
to  be  only  ca  2-3  even  in  the  most  overdamped  solvents  in  Table  I,  rather  than 
the  tenfold  or  so  predicted  from  the  dielectric-continuum  approach. 

Whither  Solvent  Friction? 

Given  this  deduction,  it  is  of  interest  to  consider  likely  factors  that  may 
account  for  the  apparent  breakdown  of  the  usual  continuum  model  of  dielectric 


17 


friction  for  activated  ET  processes.  Even  within  the  confines  of  dielectric- 

continuum  solvent  models,  there  are  good  reasons  to  doubt  the  general 

quantitative  relevance  of  to  ET  barrier-crossing  dynamics.  The  commonly 

anticipated  presence  of  faster  dissipative  solvent  relaxation  modes,  observed 

31 

directly  in  seme  TDFS  experiments,  should  accelerate  the  net  overdamped  barrier 
11  32 

crossing.  ’  This  behavior,  most  clearly  discerned  for  the  metallocene 

29  33 

reactions  in  methanol  and  other  non-Debye  media,  '  may  well  occur  more 
generally.  Furthermore,  unlike  the  solvent  inertial  motion  noted  above,  the 
longitudinal  overdamped  dynamics  associated  with  shore-range  solvation  are 
anticipated  to  differ  significantly  or  even  markedly  from  the  continuum  (zero- 
wavevector)  description  that  offers  rL_l  as  the  relevant  parameter.  This  is 
clearly  a  central  issue  since  short-range  solvation  should  inevitably  provide  the 
major  contribution  to  the  solvent  dynamics  as  well  as  the  reorganization 
energetics . 

A  number  of  recent  theoretical  and  simulational  studies  have  endeavored  to 

ascertain  the  characteristics  of  such  molecular-level  solvation  dynamics.  While 

limited  so  far  to  relatively  low-friction  media,  molecular  dynamics  simulations 

have  proved  to  be  a  particularly  insightful  approach  (refs.  39-47  are 

representative).  A  common  finding  from  these  studies  is  the  presence  of  a  rapid 

initial  solvation  response  to  the  solute  charge  change ,  approaching  that 

28 

anticipated  for  inertial  solvent  motion  followed  by  slower  dynamics  more 
characteristic  of  overdamped  relaxation.  As  pointed  out  by  Carter  and  Hynes, 
the  former  feature  is  missing  from  conventional  analytic  theories  of  solvation 
dynamics,  although  recent  treatments  have  endeavored  to  include  it.*^  Since  the 
inertial  component  can  typically  constitute  a  large  (>  50%)  fraction  of  the 
simulated  time-correlation  function,  such  rapid  molecular  dynamics  might  well  be 
anticipated  to  provide  an  important  or  even  dominant  contribution  to  ET  barrier- 
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crossing  dynamics.  Unfortunately,  however,  there  are  at  present  no  theoretical 
treatments  of  this  issue  that  yield  clearcut  numerical  predictions.  It  also 
remains  to  be  seen  if  such  rapid  short-range  solvation  dynamics  will  be 
characteristic  also  of  ostensibly  higher-friction  media  where  overdamped  motion 
is  conventionally  regarded  as  the  dominant  factor  in  the  barrier-crossing 
dynamics . 

Some  experimental  information  on  the  disparities  between  molecular 

solvation  dynamics  and  the  dielectric-continuum  predictions  has  also  been 

48 

obtained  from  suitably  ultrafast  TDFS  measurements.  The  data  indicate  the 

presence  of  relaxation  components  that  are  significantly  more  rapid  as  well  as 

slower  than  anticipated  from  the  dielectric-continuum  treatment.  The  solvents 

studied  include  anticipated  overdamped  media  such  as  dimethylsulfoxide  and 
48 

dimethylformamide .  While  limited  inevitably  by  the  pulse  laser  risetime  (ca 

0.1  ps),  these  findings  suggest  the  presence  of  unexpectedly  rapid  molecular 

(presumably  short-range)  solvent  relaxation.  Similarly  to  non-Debye  continuum 

circumstances,  such  faster  molecular-based  relaxation  components  are  anticipated 

to  dominate  the  net  barrier-crossing  frequency  for  adiabatic  activated  ET 
11,29,32 

processes. 

Consequently,  then,  the  observation  that  j/n  -  r£l  for  near-adiabatic  ET 
reactions,  as  exemplified  in  Table  I,  may  arise  from  a  fortuitous  (albeit  rough) 
compensation  between  the  influence  of  short-range  solvation  dynamical  components 
that  tend  to  be  more  rapid  than  long-range  (continuum)  solvent  relaxation, 
offset  by  the  relative  inefficiency  by  which  round-top  barriers  can  be  surmounted 
in  the  presence  of  overdamped  motion.  This  deduction  bears  a  resemblance  to  the 
conclusion  of  a  recent  MD  simulation  study  of  ET  barrier-crossing  dynamics  in  a 
model  polar  solvent,  where  the  vn  values  were  found  to  fall  only  marginally 
(twofold  or  less)  below  the  TST  prediction,  and  then  only  for  barrier-top 
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curvatures  corresponding  to  considerable  electronic  coupling  (Hu  - 
20  kJ  mol"1)  ^  (The  latter  finding  implicates  the  importance  of  frictional 
retardation  in  the  barrier-top  region,  rather  than  in  the  wells,  &  la  the  Grote- 
Hynes  description.^^)  However,  the  simulated  solvent  employed  in  ref.  46  is  only 
marginally  overdamped,  dq/tI1  -  2.5  (i.e.  comparable  to  acetonitrile  or  acetone 
in  this  respect) ,  so  that  the  relevance  of  the  findings  to  ostensibly  higher- 
friction  media  is  presently  unclear. 

Nonetheless,  in  a  more  general  vein  the  deduction  of  smaller  frictional 

retardation  effects  than  anticipated  from  conventional  zero-wavevector/Debye 

treatments  is  perhaps  unsurprising  given  the  anticipated  occurrence  of  frequency- 

dependent  solvent  friction  in  the  context  of  the  Grote-Hynes  description  of 

2b  1 1 

activated  barrier  crossings.  ’  This  treatment  predicts  that  the  barrier¬ 
crossing  dynamics  can  be  dominated  by  higher-frequency  (i.e.  faster  relaxation) 
components  of  the  molecular  solvation,  especially  when  the  barrier  top  is 
relatively  sharp.  Consequently,  the  actual  rate  retardation  in  more  overdamped 
solvents  may  be  muted  substantially  in  comparison  with  that  anticipated  on  the 
basis  of  the  major,  yet  often  slower-relaxing,  components  of  the  overall  friction 
kernel  that  are  commonly  considered  in  Debye-continuum  treatments  of  dielectric 
friction. 

In  conclusion,  one  is  left  with  an  uneasy  feeling  of  only  partial 
understanding  of  the  problem  being  achieved,  tempered  with  substantial 
uncertainty  as  to  the  underlying  physical  solvation  factors  responsible  for,  as 
well  as  the  numerical  importance  of,  solvent-friction  effects  in  activated  ET 
processes.  A  major  hurdle  to  the  exploration  of  solvent  friction  in  conventional 
(i.e.  activated)  chemical  process  clearly  remains  the  extraction  and  diagnosis 
of  dynamical  factors  in  the  inevitable  presence  of  the  other  obfuscating  effects 
that  abound  in  most  reaction  systems.  The  above  discussion  has  tried  to 


20 


emphasize  that  even  for  bimolecular  processes  useful,  at  least  semiquantitative , 
information  on  solvent-dependent  dynamical  factors  can  be  obtained  for  favorable 
cases.  Further  experimental  efforts  for  ET  processes,  however,  would  perhaps  be 
most  insightfully  directed  to  intramolecular  systems.  Such  data,  particularly 
with  information  from  HD  simulations  for  a  wider  range  of  solvating  media,  may 
well  spawn  a  new  chapter  in  our  molecular-level  understanding  of  solvent-friction 
effects  in  activated  chemical  processes. 
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TABLE  I  Comparison  Between  Solvent-Dependent  Barrier-Crossing 

Frequencies  for  Selected  Metallocene  Self  Exchanges  and  Predicted 
Inertial  and  Overdamped  Solvent  Frequencies 


*.1^  s 


SOLVENT  B 

w0/2ir  b 
s'1 

s  1 

uB  d 

S^1 

Cp2Co+/0  f 

Cp*Co+/0  8 

Cp2Fe+/0  h 

Acetonitrile 

2xl012 

4xl012 

1 . 2xl012 

5xl012 

2.5xl012 

3xl0u 

Acetone 

l.OxlO12 

3 . 5xl012 

6.5xlOu 

3xl012 

2xl012 

4x1 01 3 

Water 

7xl012 

1 . 9xl012 

7.5xl0n 

3 . 3xl012 

Nitrome thane 

1 . 6xl012 

4.5xl012 

8xl0n 

3xl012 

4xl0n 

DMSO 

1 . 5xl012 

5xl0n 

1.7xl0u 

l.OxlO12 

1.5xlOu 

Benzonitrile 

6 . 5xl0n 

2xl0u 

5 . 5xl010 

7xl0n 

3xl0n 

2 . 5xlOu 

Nitrobenzene 

-  6xl0n 

2xl0n 

6.5xl010 

1.5xlOu 

2.5xlOn 

TMU 

8xl0n 

1 . 7xl0u 

6xl010 

3xlOn 

HMPA 

-  6xl0n 

l.lxlO11 

4xl010 

2 .  OxlO11 

PC 

1 . 5x10 12 

(4xlOu) 

2xlOn 

2 . 5xl012 

l.SxlO12 

3 . 5xlOn 

Methanol 

1 . 7xl012 

(1.3xl0u) 

4xlOn 

1 . 5xx012 

4.5xlOn 

DMSO  -  dimethylsulfoxide ,  TMU  -  tetramethylurea ,  HMPA  —  hexamethylphosphoramide ,  PC 
-  propylene  carbonate. 

b  Predicted  adiabatic  TST  barrier-crossing  frequency  in  given  solvent,  from  Eq.  (3). 
(See  footnotes  to  Table  II  and  Fig.  4  of  ref.  29  for  calculational  details.) 

Inverse  longitudinal  relaxation  time  of  pure  solvent,  as  extracted  from  dielectric 
loss  spectra  (from  Table  I  of  ref.  22).  Values  given  in  parentheses  refer  to  major 
relaxation  for  solvents  also  exhibiting  higher-  frequency  dispersion  (non-Debye 
behavior) . 

d  Adiabatic  barrier-crossing  frequency,  calculated  from  Eq.  (10)  by  using  dielectric 
continuum  model,  assuming  H12  -  2 . 1  kJ  mol-1  (see  text  and  ref.  29  for  details). 

£ 

Barrier-crossing  frequencies  for  metallocene  redox  couple  in  given  solvent,  as 
extracted  from  rate  constants  and  barrier  heights  (obtained  from  nmr  line  broadening 
and  near-infrared  optical  measurements,  respectively)  by  using  Eqs .  (2a)  and  (2b). 
See  Table  I  of  ref.  22  for  experimental  data,  and  ref.  21  for  procedural  details. 
Value  of  precursor  stability  constant,  Kp,  taken  as  j.25  M-1.  Barrier  heighcs  for 
Cp2Co+/0  and  Cp|Co+/0  reactions  diminished  by  2 . 1  and  3 . 1  kJ  mol'1,  respectively,  to 
allow  for  electronic  coupling  (ref.  22). 

^  Cp2Co+/0  -  decamethyl(cobaltocenium  -  cobaltocene) . 

8  Cp2Co+/0  -  carboxymethyl(cobaltocenium  -  cobaltocene). 

b  Cp2Fe+/0  -  ferrocenium  -  ferrocene. 


